Physical origin of the buckling in CuC^: Electron-phonon coupling and Raman spectra 
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It is shown theoretically that the buckling of the Cu02 planes in certain cuprate systems can 
be explained in terms of an electric field across the planes which originates from different valences 
of atoms above and below the plane. This field results also in a strong coupling of the Raman- 
active out-of-phase vibration of the oxygen atoms (Bi g mode) to the electronic charge transfer 
between the two oxygens in the Cu02 plane. Consequently, the electric field can be deduced from 
the Fano-type line shape of the B\ g phonon. Using the electric field estimated from the electron- 
phonon coupling the amplitude of the buckling is calculated and found to be in good agreement 
with the structural data. Direct experimental support for the idea proposed is obtained in studies of 
YBa2Cu306+a; and P^S^Cai-^Y^C^Os with different oxygen and yttrium doping, respectively, 
including antiferromagnetic samples. In the latter compound, symmetry breaking by replacing Ca 
partially by Y leads to an enhancement of the electron-phonon coupling by an order of magnitude. 

PACS numbers: 74.72.-h, 63.20.Kr, 78.20.Bh, 71.I0.-w 



I. INTRODUCTION 

The physical origin of the buckling of the Cu02 plane in certain high temperature superconductors, where the 
oxygen atoms are placed outside of the plane of the copper atoms, is a long standing problem. There are two possible 
explanations: (i) the surroundings of the CuC>2 plane is not symmetrical, thus an electric field perpendicular to 
the plane acting on atoms with different charges deforms that plane, and (ii) the energy associated with chemical 
bonds (especially the bond between the oxygen 7r-orbitals perpendicular to the plane and the coppers) results in 
a spontaneous symmetry breaking. The present paper provides direct experimental evidence to support the first 
possibility in accordance with the theory. 

R has been suggested previously |lj that in the presence of a perpendicular electric field a charge transfer between 
the two different oxygens 0(2) and 0(3) invokes a deformation of the type of the B\ g Raman-active phonon at 
330cm -1 , where the two oxygens move in an out of phase manner perpendicular to the plane (Fig. la). The resulting 
electron-phonon coupling was then used to interpret the Fano interference in the Raman spectra of YP^CuaOg+z 
(Y-123), where the light is scattered by both the B\ g phonon and the electrons involved in the charge transfer in the 
plane described above. For a Cu02 plane in a more symmetrical environment (where the electric field is negligible) 
such a linear coupling does not exist since the energy correction due to the deformation of the chemical bonds is at 
least quadratic in the displacement of the oxygens. Hence, only much weaker two-phonon processes contribute. 

The coupling strength and, consequently, the value of the electric field perpendicular to the plane can be determined 
from the Fano line shape of the B\ g phonon using the expressions derived in Ref. Q). In order to follow closely the 
doping dependence of the line shape we will use here, in contrast to [[lj, a more realistic three band model for the 
electronic system. It will be described in some detail thus augmenting an earlier publication which gave just a 
brief summary [0. The Raman- active continuum which interferes with the B\ g phonon is assumed to originate 
from strongly relaxing electrons in a single conduction band such as proposed in Refs. |^-|5| and will be introduced 
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FIG. 1. The unit cell of the Cu02 plane in Y-123 shown with the atomic displacements corresponding to the Bi g (a) and 
Ai g (b) phonons. The electric field E perpendicular to the planes is due to the asymmetric environment (Ba 2+ above, Y 3+ 
below), and causes a static deformation of the A lg type, usually referred to as buckling (also called dimpling as in Ref. ]H|), 
8p and the double arrows denote the charge transfer which accompanies the lattice vibration. In Bi-2212, Ba 2+ is replaced by 
Sr 2+ and, at optimal doping, Y 3+ by Ca 2+ . Then the electrical field E is negligible. 



phenomcnologically as in y. Recently, Sher man and coworkers considered the coupling of the B\ g phonon to non- 
interacting Bloch electrons including intra and inter-band intermediate states in order to understand the phonon 
intensity Here, the focus is placed on the interference of the phonon and the electronic continuum which are 
excited by the light through two initially independent channels. In our approach the direct coupling of the phonon 
to the light is parameterized and not explicitly considered. Our goal is to better understand the interaction between 
the phonon and the conduction electrons. Therefore, both Y-123, where the surroundings of the Cu02 plane is highly 
asymmetric, and I^S^CaC^Og (Bi-2212), where the plane is placed in a more symmetrical position with respect 
to the charges of the nearby atoms, will be studied. 

The B\g phonon in metallic cuprate systems has actually been studied intensively by several groups The 
main emphasis was placed on the anomalies at the transition to the superconducting state. It was shown, for instance, 
that the change in frequency of the B\ g phonon at T c vanishes already at moderate underdoping || pointing at a 
considerable change in the elect ron-phonon coupling. While the phonon renormalization at T c has been the subject 
of numerous theoretical treatments |10[ , no attempt has been made yet to link the doping dependence of the coupling 
to a realistic band structure, i.e. a microscopic model. For a comparison with the model predictions, a systematic 
experimental study of the phonon covering the complete doping range from antifcrromagnetic (AF) insulating to 
slightly overdoping in comparable, high-quality samples is essential. Therefore, we present not only data on differently 
doped metallic samples but also new data on undoped antiferromagnetic ones in order to get an idea of the intrinsic 
line width and to study in detail the influence of a reduction of the continuum scattering on the line shape of the 
phonon. 

As already shown earlier j^,|?J the B\ g phonon has a reasonably Lorentzian line shape for the symmetric compound 
Bi-2212 indicating the lack of any substantial coupling between that phonon and the charge transfer between the two 
oxygens in accordance with the absence of an electric field. This striking difference to the data of Y-123 || can be 
considered as strong evidence in favor of the linear (involving one phonon) electron-phonon coupling due to the electric 
field. In order to complete this argument, we gradually replace Ca 2+ by Y 3+ in the Bi-2212 compound where the 
doping most likely breaks locally the approximate reflection symmetry through the Cu02 plane, thereby making this 
material similar to Y-123. In fact, the phonon becomes more asymmetric and the electron-phonon coupling increases 
by an order of magnitude However, this increase does not at all result in an enhancement of the superconducting 
transition temperature: T c is rather reduced on substituting Ca with Y and follows the usual T c vs. doping curve. 

Finally, the value for the electron-phonon coupling obtained from the analysis of the line shape is found to be 
in excellent agreement with the value calculated from the magnitude of the electric field due to the experimentally 
observed buckling |l^] and the restoring force calculated from the frequency of the A\ g phonon (Fig. lb). The 
coupling depends on changes in the doping via the proximity of the van Hove singularity in the density of states to 
the Fermi level Ep. 

The paper is organized as follows: In section II the theoretical model is described in some detail and the relevant ex- 
pressions are derived. In section III and IV experimental details and results are being given. Finally, the consequences 
are discussed. 
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II. THEORETICAL MODEL 



Much of the development of the crystal field model was given in Ref. pj] , and details of the calculations are found 
therein. We start with this and only point out where changes are made. The model constructed in Ref. (I) started from 
a three band model where only the Cu-0 hopping was taken into consideration. However, one must add 0-0 hopping 
in order to represent the curvature and centering of the Fermi surface as revealed from angle-resolved photo-emission 
experiments. Therefore we augment the results from jl| to include direct 0-0 hopping. We start by first considering 
a simple 3 band model for the Cu02 plane with Cu-0 hopping amplitude t and 0-0 hopping amplitude t': 

H a = eJ2 *C V* + t £ [P S bl„a n ^8 + h.c] + f£ £ ^^^V' W 

a^,S n ' a (6,6') 

where a creates an electron with spin a at a, copper lattice site n, while a n a § annihilates an electron at one of 
the neighboring oxygen sites n + 6/2 determined by the unit vector 6 assuming the four values, (±1,0) and (0, ±1). 
An oxygen atom between the two copper atoms at sites n and n + 6 is labeled by either (n, 6) or (n + 6, —6). As 
in jl| , e — Ed — E p is the difference of the Cu and O site energies and Pg = ±1 depending on whether the orbitals 
(with real wavefunctions) have the same or opposite sign at the overlap region. Assuming Cu d x i_ y i and O p orbitals 
P ' g = —P§, and we can choose P(i,o) = 1 an( i P(o,l) = — 1- Lastly, Pg g> denotes the overlap sign between an oxygen 

orbital at site n + 6/2 with a neighboring oxygen orbital at site n + 6 /2. By our above convention these overlaps 
take the values P^ y = P'_ x _ y = 1, P x _ y = P- xy = — 1, respectively. After Fourier transforming, the Hamiltonian 
now reads H° — J2 kcr H^ a , where 

H k,cr = £ &k, CT &k,<T + {^LcrKi^aMk) - a 2/,k,<A(k)] + h.C.} + 4[<4,k,<7 a 2/,k,<7 + h.C.], (2) 

with the prefactors 

i a (k) = 2isin(afc Q /2), i'(k) = -At' sin(afc x /2) sin(ak y /2). (3) 
We can then easily diagonalize Eq. (Q) and obtain the 3 bands: 

< CT = ^^(k)4 k)CT ^, k>CT , (4) 



where (3 assumes the values +,— , and for the anti-bonding, bonding, and non-bonding bands, respectively. The 
energies are then given by solving a third order root equation and are 

E+(k) = s+ (k) +s_(k) +e/3 (5) 

E-(k) = -\{s + (k) + S _(k)] + s/3 + ^> + (k) - Mk)] 

E Q (k) = -\[s + (k) + s _(k)] + e/3 - ^[ S+ (k) - s _(k)] 

with 

s± (k) = (r(k) ± v/<7 3 (k)+r 2 (k)) i (6) 

g(k)=-i[^(k)+^(k)+t«(k)]- e a /9 

r(k) = £ -[tl(k)+t 2 v (k) - 2t' 2 (k)} -t'(k)t x (k)t y (k) +e 3 /27 

Since g 3 (k) + r 2 (k) is negative for all k, s!J_(k) = s_(k) and the energies are of course real. 

The parameter values used, e = leV, t — 1.6eV, and t'/t — 0.35, are similar to those used in Ref. [jll]]. The 
dispersion for the three bands is plotted for symmetry k points along the Brillouin zone (BZ) in Fig. || using these 
parameters and choosing the chemical potential such that the filling (n) — 0.85 for both spins. The upper (middle, 
lower) curve represents the band dispersion for the anti-bonding (non-bonding, bonding) band, respectively. Setting 
the direct 0-0 hopping integral t' = would yield a dispersionless non-bonding band whose energy equaled the 
oxygen site energy, while the other bands recover the dispersion calculated in Ref. M . The inset of Fig. @ shows the 
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FIG. 2. Band dispersion for the 3-band model with nearest and next - nearest - neighbor hopping for the parameters listed 
in the text corresponding to optimally doped YBa2Cu3C>7. The upper (middle, lower) curve corresponds to the anti-bonding 
(non-bonding, bonding) band, respectively. Inset shows the doping dependence for the upper band only in the vicinity of the 
van Hove point near (n, 0). The upper (middle,lower) curve represents over- (optimally, under-)doped systems, respectively. 



dispersion of the anti-bonding band in the vicinity of the (it, 0) point in the BZ evaluated for different fillings which 
we label as over-doped ((n) = 0.8), optimally-doped ((n) = 0.85), and under-doped ((n) = 0.875), corresponding to 
the upper, middle, and lower curves in the inset, respectively. The filling changes the distance the flat part of the 
dispersion (van Hove point) is away from the Fermi level Ef, thus changing the density of states at Ef- Through 
this change in the density of states we attempt to understand the doping dependence of the coupling between the B\ g 
optical phonon mode and the conduction electrons. 

At this point we can carry over the results from Ref. |l| with little loss of generality. As in Ref. [jlj we consider only 
a reduced one band model appropriate for near half filling and take only the upper band E + (k) into account. The 
electron - phonon reduced Hamiltonian for the J5 lg q = phonon mode was given as 



H el . ph ^ 



^^ 5 (k)4^ k>CT [c + C t] 5 (7) 



where N is the number of sites, dk,(r annihilates an electron of spin a and momentum k, and creates a B\ g phonon 
mode of wavevector q = 0. The coupling constant <?(k) of the B\ g mode to an electron with momentum k was 
evaluated in Ref. JlJ and is given by 



S(k) = eE \j^^^\ I 2 - I I 2 ]. (8) 

where M is the oxygen mass, w_b 1s is the phonon frequency, and E is the z- component of the electric crystal field 
which for simplicity we take to be equal at both the 0(2) and 0(3) sites. The functions (j> xy are the amplitudes of 
the oxygen wave functions in the tight-binding wave functions of the upper band which arise from diagonalizing the 
three band model. In the absence of 0-0 hopping, these functions were given in Here, 0-0 hopping is included 
yielding 

= J^[E + (k)t X (k) - t'(k)ty(k)], 

<Mk) - j^r[E + (k)t y (k) - t'(k%(k)], 

with iV 2 (k) = [E 2 + (k) - t' 2 (k)} 2 + [E+(k)t x (k) - ^(k)t'(k)] 2 + [E+(k)t y (k) - t x (k)t'(k)} 2 . (9) 

We now are in a position to calculate the full Raman response. Again we can utilize the results derived in Ref. 
. The generalized form of the Breit-Wigner or Fano lineshape describes the scattering of the light by the coupled 
electron-phonon system. The (independent) coupling strengths of the light to the phonon and to the electronic 
continuum with symmetry A will be denoted by g p - Pi \ and j\, respectively. The bare phonon when decoupled from 



the conduction electrons is characterized by the energy lu\ and the intrinsic damping due to e.g. anharmonic phonon- 
phonon interaction by T\. The electronic continuum in symmetry channel A will be described by the electronic 
susceptibility \\- It is clear that it cannot originate from a non interacting gas formed by Bloch electrons; it rather 
should be traced back to strong renormalization effects such as electron-electron interaction || , spin-fermion || or to 
scattering on static impurities Q . A parameterized form for \\ is used which has been derived by one of us when 
electron dynamics is considered in materials with nesting properties along the Fermi surface p]: 



LOT, 



(10) 



where 



A = A 



Cu = Lom*{u>)/m, 

with m*(id)/m = 1 H In 

7T 



(11) 



Here Np is the density of states for both spins at the Fermi level, m, m* is the bare, renormalized electron mass, 
respectively, = 1/r — 1 /t\ is the channel-dependent impurity scattering rate reduced by vertex corrections Q] 
and r is the lifetime of the electrons, a, 0' and lo c are constants determined by a fit to the electronic continuum. 
The cut-off w c is essentially the band width, /3' measures the interaction between the electrons and can vary only 
between 3.3 and 4.2. They both will be fixed as indicated in Table I. a < 1 as long as f3' < 3.7 ||. In addition, 
a ■ f3' determines the temperature dependence of the normal-state d.c. resistivity and should therefore be close to 1.5. 
This means that there exist strong theoretical constraints for metallic samples. 1/t^ is a free parameter to get an 
optimal representation of the continuum. This form for the Raman susceptibility provides an adequate description to 
the continuum in the normal state of both Y-123 and Bi-2212 (see Ref. As the continuum depends only weakly 
on frequency it turns out that the choice of the parameters mildly influences the Fano line shape. The full Raman 
response measured in channel A for the interacting electron-phonon system reads JlJ : 



Wa 



LU + LU a 



7aXaH 



[1 + AH//3] 2 



WA 



UJ + UJ a 



(12) 



Here, the first factor on the right hand side is the Lorentzian line shape of the phonon modified by the Fano interference 
and the second one describes the interplay between the scattering of light on phonons and electrons. The details depend 
on the relation between the coupling constants g p - p \, which has been studied by Sherman and coworkers ||, 7a, and 
g\. The latter one describes the interaction of the phonon with the conduction electrons and can also be calculated 
directly from the band structure and the electric field by averaging the coupling g(k) given in Eq. (^J) over the Fermi 
surface. For B\ g symmetry we obtain 



-E 



9(k) 



de(k) ' 



(13) 



where / is the Fermi function and the factor of two accounts for both spin directions fli3| . For convenience the 
dimensionless quantities \{u>) and (3 were introduced in Eq. (O) which depend on the coupling constants as follows 



and 



2g P - P ,B lg gB lg 



lB lg UB lg 



(14) 



(15) 



X(u>B lg ) is the usual electron-phonon coupling constant. Here the explicit dependence of A(a>) on to through \' Bi (ui) 
has been dropped as the real part of the electronic response x'b 19 ( w ) m the normal state turns out to be a smooth 
function of ui close to tOB lg - Using these expressions we obtain the resonance frequency of the phonon Cjb 1cj 
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£! 19 =< 9 [1-A(wb 1s )] (16) 

which depends only on the electron-phonon coupling. The position of the antircsonance io\ can be given as 

w2=w| la [l + /J]. (17) 

Through (3 the antiresonance frequency depends also on ratio of the photon-phonon coupling g p ~ p ,B lg to the projected 
photon-electron coupling jg 1 . The magnitude of these parameters is not important and enters as a multiplicative scale 
factor determining the absolute cross section in the same way as the density of states Np- Finally, electron-phonon 
coupling leads to a larger line width 

r Bl >iO = IV, + g 2 Blg XB lg (uB lg )/N F , (is) 

neglecting strong anisotropy of the electron-electron or electron-impurity self energies. 



III. SAMPLES AND EXPERIMENTAL 

The Y-123 crystals were grown in BaZrC>3 (m) which has been shown to be completely inert and to facilitate the 
preparation of samples with a purity of better than 99.995%. All crystals were postannealed in pure oxygen and 
quenched. Temperatures and oxygen partial pressures were adjusted according to the calibration of Lindemer et al. 
p5| . The resulting oxygen concentrations were approximately 6.0 6.5, 6.93, and very close to 7.0 for the samples we 
call antiferromagnetic (AF) insulating, under-doped, optimally doped and over-doped, respectively, in the following. 
The magnetically determined respective T c values and transition widths of the superconducting samples were 53.5 A" 
(AT C = 3 A), 91.5 (AT C =0.3 A), and 87.0 (AT C = 1.0 A). 

The Bi-based samples, Bi2Sr2(Cai_ 2; Y 2 ;)Cu208 (Bi-2212), were prepared in ZrO crucibles. In optimally doped 
crystals (without Y) the resistively measured T c was generally above 91 A with AT C < 2A. If Ca 2+ is replaced by 
Y 3+ (see Fig. |l|) holes in the CuC>2 planes are filled in and T c is reduced. At the same time the environment of 
the planes becomes electrically asymmetric. The samples we used contained 38% Y and 100% Y, respectively. The 
superconducting sample (38%) was well in the underdoped range of the phase diagram with a T c of 57 A (AT C = 5 
A). Bi-2212 is not stable at the stoichiometric composition and there is always excess Bi in the crystals (typically 
Bi 2 .i instead of Bi 2 ) fT|Jl7| which is found predominantly in the Ca (Y) position. In addition, a small amount of Sr 
may be replaced by Ca or Y. Therefore, the environment of the Cu-0 plane is by far less ordered than Y-123. In 
general, it cannot be expected for Bi-2212 that a crystal quality comparable to the one of Y-123 is obtainable. By 
partially substituting Bi with Pb the modulation along the crystallographic b-axis can be changed or even suppressed 
completely fl9| . In the sample we used the Laue pattern still showed an indication of a twofold symmetry but it is 
reasonable to assume that the distortion is smaller than in Pb-free samples. 

The experiments were performed in back-scattering geometry using a double monochromator with single- channel 
detection and the resolution set at 8 cm -1 . This is sufficient as the total widths of the narrowest lines are roughly 
twice as large and add quadratically to the resolution. For excitation the Ar + line at 476 nm was selected. The 
maximal power was 4 mW in order to keep the laser- induced heating below 15 A. The beam was focused to a spot of 
approximately 50 x 150/im 2 . For the study of the excitations we were interested in the polarizations of the incoming 
and outgoing photons were always parallel to the planes. The coordinate system is locked to the Cu-0 bonds with x 
= [100], x' — [110], etc. All symmetries refer to a tetragonal pointgroup. B\ g phonon and continuum are projected 
out with x'y' polarization. 



IV. RESULTS 

Results for Y-123 obtained at B\ g symmetry {x'y') are plotted in Fig. |^. All spectra are divided by the Bose- 
Einstein thermal function in order to get the response Imx as given in Eq. (12). As a result of doping the shape of 
the B\g phonon at approximately 330 cm -1 changes considerably. For the AF and slightly doped sample it is narrow 
and close to a Lorentzian. When carriers are added the line broadens and becomes more asymmetric exhibiting 
a Fano-type dependence on frequency. The oscillator strength of the phonon decreases while the B\ g continuum 
gains intensity pOjj . As a consequence the intensity ratio Iphonon to I C ontinuum decreases by a factor of at least 20. 
We emphasize that the low energy continuum vanishes almost completely in the AF crystal, hence supporting the 
interpretation in terms of charge excitations. 

In optimally doped Bi-2212 (Fig. [| (c)) the B\ g phonon line is much weaker while the continuum has roughly the 
same cross section as comparable Y-123 (Fig. || (c)). The bigger line width comes at least in part from inhomogeneous 
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FIG. 3. Big spectra for differently doped YBa2Cu306+^ as indicated. Shown are experimental data (after division by the 
Bose factor) and theoretical results from the crystal field model (Eqs. (Jl^)-(Jl^)). The parameters used are defined in the text 
and compiled in Table I. 
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FIG. 4. Big spectra for differently doped (Bi,Pb)Sr2(Ca,Y)Cu208 as indicated. Shown are experimental data (after division 
by the Bose factor) and theoretical results from the crystal field model (Eqs. (|l^)-(|l8|)). The parameters used are defined in 
the text and compiled in Table I. 
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TABLE I. Summary of fitting parameters used in Figs. |3| and |4j. T c is the superconducting critical temperature, T the 
measuring temperature. AF indicates that the sample is antiferromagnetic and insulating. The seven parameters in the third 
column are the fitted parameters. The actual Fano profile determines the parameters i^B lg , i^a, and A which depend only 
weakly on a, t% x _1 , /3', and uj c (see text). (2>B lg is calculated from uiB lg and A. The last column are fixed parameters. A scale 
factor was multiplied to Eq. ([l2]) to account for the overall magnitude of the cross section. The scale factor is proportional to 
the density of states at the Fermi energy Nf and is, as expected, very small in the AF samples. 

broadening as explained above which manifests itself in the large intrinsic line width (see Table I). If Y is doped in 
for Ca a new line at approximately 315 cm -1 starts to gain intensity very rapidly (Fig. |] (a) and (b)). A reminder of 
the line found in Y-free crystals is still seen as a shoulder on the left-hand side of the new one (Fig. [| (b)). This tells 
us that the Y is not completely homogeneously distributed in the underdoped sample. Most likely, the Y replaces 
big regions of the Ca planes, just like in Y-123 where upon oxygen doping clusters of longer and shorter fractions 
the chains build up depending on the conditions of preparation Jis] . So, in addition to changing the charge balance 
around the Cu02 plane, defects are introduced in its vicinity which are seen in both the phononic and the electronic 
properties. Apparently the intrinsic lincwidth in the AF sample is somewhat smaller than in the metallic ones. This 
comes most likely from the additional doping with Pb which is used to suppress the superstructure of the Bi02 layers 
and the related orthorhombic distortion of the unit cell. The integrated cross section is now as large as in AF Y-123 
[compare Fig. ^ (a) and Fig. ^ (a)]. This is further independent evidence for the model proposed. 

V. DISCUSSION AND CONNECTION TO BUCKLING 

Excellent fits to the data on optimally-doped YBa2Cu307_a in both the superconducting and normal states were 
given in [Q, which when interpreted supported evidence for c? a; 2_ l) 2-pairing in this material. We have extended these 
results by examining normal state data of different samples and dopings. We note that a completely new set of Y- 
123 single crystals with superior quality (see section III) is used here. 

We find that once again very good fits to the data can be obtained on four and three differently doped samples of 
Y-123 and Bi-2212, respectively. For all fits given in Figs. | and |, the values used to calculate the susceptibilities of 
nested Fermi liquid theory and the phonon lineshapes are summarized in Table I. 

Eight parameters (a, tJj , T l Bi , wg ls , ui a , A, u)B lg , and an overall scale factor) are used to provide the fits in Figs. 
H and ||. However, we do not have eight separate degrees of freedom, a, t Bi ^ , and the overall scale factor are chosen 
to reproduce the continuum for frequencies away from the phonon. t Bi determines the low frequency part of the 
continuum while a determines the high frequency part. The uncertainty in these parameters is less than 10%. In 
addition, the parameters a and (3' are in the narrow range imposed by the nFl model. || Then tliB lg is chosen to 
locate the center of the phonon line, which can be accurately placed to within a fraction of a wavenumber (i.e. much 
better than the experimental resolution) . Next A is chosen to reproduce the local asymmetry around the phonon line 
(Fano). This automatically fixes ujB lg via Eq. ( |l6| ) and sets the damping from the coupled electron-phonon system. 
Then, the intrinsic line width T B is set to account for any remaining damping needed to account for the phonon 
linewidth. An indication of the stability is the almost constant intrinsic line width in Y- 123. Indeed these samples are 
comparable and the unit cell changes only slightly. On the other hand, it is not surprising that the intrinsic line width 
is not constant in Bi-2212 since the environment of the planes and even the total crystallographic structure changes 







considerably on doping. Finally, u> a is then chosen to set the overall intensity of the phonon relative to the electronic 
continuum and to pinpoint the minimum of the spectrum on the high frequency side of the phonon. In summary, 
each parameter is independently chosen and the accuracy of the fit can be determined by the agreement of the theory 
to the experimental background and the phonon's position, intensity, asymmetry, linewidth, and local high frequency 
minimum. Therefore we are confident that the parameters used are accurate to within the 10% level. Should further 
constraints be imposed (i.e. by a more accurate theory for the continuum or the phonon-photon coupling constant), 
it would be straightforward to adapt the theory accordingly. 

We now can make a comparison to the crystal field model predictions for the electron-phonon coupling constant. 
We see from the fits that for the Y-123 compounds the electron-phonon coupling increases up to a factor of two with 
increasing doping. This can be understood in part by an increase in the density of states at the Fermi level with 
increased oxygen doping as the van Hove contribution moves closer to the Fermi level (see Fig. ||). Our calculations 
using Eq. ( Ilif ) yield the following values for A as a function of doping: for A = 0.0321 for < n >= 0.875 (underdoped, 
Fig. 3b), A = 0.0382 for < n >= 0.85 (optimally doped, Fig. 3c), and A = 0.0482 for < n >= 0.8 (overdoped, Fig. 
3d). Here we have taken T = 100K,ujB lg = 348cm _1 , M = 16mp, where rap is the proton mass, and used a value of 
the electric field which is 1.3V/A |23[ . We observe that these values of A are quite close to the values used to obtain 
the fits in Figs. || (b)-(d). Given the values of t'/t,t, and e are taken to be independent of doping, the agreement is 
quite good. The agreement could be refined once parameter choices for various levels of doping are determined via, 
e.g., fitting to experimentally observed Fermi surfaces. 

We can connect the value of the electric field used to account for the observed buckling of the Cu02 plane in Y-123. 
The electric field E perpendicular to the CuC>2 plane results in an Ai s -type static distortion of the plane since the 
charges of oxygen and copper are different. It is very hard to perform a complete calculation as the electric field 
is different at the Cu and the O sites. Furthermore all the lattice forces have to be known. For an estimation of 
the buckling, however, a simplified model is sufficient |l2j where Cu is pinned rigidly to the elementary cell and the 
oxygen moves in a harmonic potential being characterized by the frequency of the A\ g phonon at u>A lg — 435 cm -1 
pjj . The restoring force at the buckling amplitude Az must balance the electric force acting on the oxygen with 
charge q = -1.75e. Thus, qE = Mu Ai Az must hold. With the experimental value Az = 0.241 @, E = 1.53V/ A 

is obtained which is close both to the 1.3V/ A used for estimating the electron-phonon coupling strength and to 
the theoretically calculated number p3}| . On the other hand, in the case of Cu02 planes in a more symmetrical 
environment as in Bi-2212 and the infinite-layer compound CaCu02 the buckling, if it exists at all, is at least an order 
of magnitude smaller as found in structural studies p4] , p5| . 

At the same time, the electron-phonon coupling for optimally doped Bi-2212 (Fig. 4c) obtained from the fit is more 
than one order of magnitude smaller than for any of the metallic Y-123 compounds and, as a result, the weight of the 
phonon is also reduced by at least the same amount. Indeed, in our model we would expect that the electron-phonon 
coupling is significantly smaller since the crystal electric field must be much weaker in this compound than in Y-123 
(see Fig. [I]). We attribute then the small electron-phonon coupling to a much smaller local asymmetry. A weaker 
spontaneous symmetry breaking is also possible. 

To test this idea we presented measurements of a sample of Bi-2212 which has been doped with small amounts 
of Y in place of Ca (Fig. 4b). Since here the valence of Y (+3) is different from that of Ca (+2), once again the 
mirror plane symmetry is broken and we would expect a much larger electron-phonon coupling than in the undoped 
compound. The fit to the B lg data using Eq. ( [l2|) is given in Fig. ||(b) and yields A = 0.0124. The coupling constant 
is indeed increased over that of the optimally doped compound by an order of magnitude when Y is introduced and 
in accordance the intensity of the phonon line is also essentially enhanced. This trend continues in the AF sample 
Fig. 4(a), and the phonon cross section is now similar to that in Y-123. 

The intensity of the continuum in Bi-2212 Y does not approach zero as one would expect for an insulator. In- Y-123, 
too, the continuum is not vanish identically but is significantly smaller. At present, we do not know the origin of this 
intensity. It can be either photoexcited carriers as observed in semiconductors p2[ , an effect of remaining carriers, 
or coupling of the phonon to the spin system p6[ . As a consequence the coupling A is not as small, though reduced 
considerably, as to be expected and actually observed in Y-123. This again gives strong support for the crystal field 
coupling model as the driving source of electron-phonon coupling. 

VI. SUMMARY AND CONCLUSIONS 

In summary, studying Y-123 with different doping levels we have shown that the breaking of the local reflection 
symmetry through the Cu02 planes due to the position and charges of the neighboring atoms results in an electric 
field across the planes which is sufficiently strong to produce both the observed buckling and the strong electron- 
phonon coupling being responsible for the Fano-type line shape of the respective B\ g mode. In order to check this 
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idea experiments were performed on Bi-2212 with and without Y doping. While the sample without Y shows very 
weak electron-phonon coupling, the interaction is enhanced by an order of magnitude and becomes comparable to 
the one in Y-123 if the local reflection symmetry is broken by replacing part of the Ca by Y. (For Bi-2212 and 
CaCu02 there is no experimental indication for spontaneous symmetry breaking.) As doping Bi-2212 by Y results in 
a change of T c from T c = 91.5K to T c = 57K accompanied by a large increase of the coupling A, the u> = 330cm -1 
phonon can not play an important role for the superconductivity, in agreement with the conclusion of Savrasov and 
Andersen [ p7| . We conclude that the symmetry breaking due to the composition can explain both the buckling of the 
planes and the strong linear electron-phonon coupling in the two types of material. One should keep in mind that 
spontaneous symmetry breaking can occur also in other materials like La2_a;Sr x Cu04 and determine the properties. 
Further experimental and theoretical studies of CaCu02 would be useful. It is of interest in the present context that 
the so called infinite layer material CaCu02 does not show any sign of buckling , consistently with the symmetric 
position of the Cu02 plane. 
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